Electric field effect modulation of transition temperature, mobile carrier density and 
in-plane penetration depth in NdBa2Cu307 ^ thin films 
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We explore the relationship between the critical temperature, Tc, the mobile areal carrier density, 
n2D, and the zero temperature magnetic in-plane penetration depth, Aab(O), in very thin underdoped 
NdBa2Cu307_i films near the superconductor to insulator transition using the electric field efi'ect 
technique. Having established consistency with a Kosterlitz-Thouless transition we observe that 
Tkt depends linearly on n2D , the signature of a quantum superconductor to insulator transition in 
two dimensions with zV = 1, where z is the dynamic and V the critical exponent of the in-plane 
correlation length. 

PACS numbers: 71.30.-|-h, 74.78.Bz, 74.72.-h, 74.75. Ha 

The electronic properties of high Tc superconductors 
are critically determined by the density of mobile holes 
as illustrated by their generic phase diagram in the tem- 
perature dopant concentration plane. Understanding the 
physics of this phase diagram, in particular at the two 
T = edges of the superconducting dome, has emerged 
as one of the critical questions in the field of cuprate 
superconductors. In the underdoped regime, where a su- 
perconductor to insulator transition occurs, there is an 
empirical linear relation (the so-called Uemura relation) 
between Tc and A~j^(0) 1], where Xab is the in- plane Lon- 
don penetration depth. If, in the underdoped limit, a 
two dimensional (2D) quantum superconductor to insu- 
lator (QSI) transition occurs, such a relation between Tc 
and 1/A^(0) is expected This however necessitates 
the occurrence of a 3D to 2D crossover as the under- 
doped limit is approached with a diverging anisotropy 
7 = Xc/Xab- In the case of a 3D QSI-transition, one 
would expect Tc oc Xa,b,c where z is the 

dynamic critical exponent of the transition Recent 
penetration depth measurements on YBa2Cu306+a; sin- 
gle crystals (where 7 apparently saturates at low 
doping levels) and films suggest a 3D-transition with 
22:/(z-fl)~l|^,a result different from the Uemura re- 
lation. Furthermore, an empirical relation involving the 
normal state conductivity and extending up to optimum 
doping was recently proposed by Homes et al. [7|. How- 
ever, all these studies on the relationship between Tc and 
Xa,b,c (0) stem from samples where the doping level was 
modified by chemical substitution Q,|lllElEllSli- 

In this letter we use the electrostatic field effect to 
tune the mobile areal carrier density, n2D) in ultrathin, 
strongly underdoped NdBa2Cu307_5 (NBCO) films and 



explore the intrinsic relation between Tc, n^o and A~j^(0) 
close to the superconductor to insulator transition. The 
electrostatic field effect technique, which allows the car- 
rier density in a given sample to be changed without af- 



fecting the chemical composition, thus avoiding sample 
to sample variations and substitution induced disorder, 
appears to be an ideal method to elucidate the intrinsic 
relationships between the key physical parameters Q. 

In a classical oxide field effect configuration, an electric 
field is applied across a gate dielectric in a heterostructure 
made of an oxide channel and a dielectric layer. The 
density of carriers is modified at the interface between the 
oxide channel and the dielectric, changing the electronic 
properties of the oxide channel. This technique allows us 
to induce in thin oxide superconducting films rather large 
Tc modulations In the 3-4 unit cell thick films used 
in this study, Tc modulations of more than 10 K have 
been achieved, corresponding to induced charge densities 
of the order of 0.7 • lO^** charges/cm^. 

The field effect device used in this study, described 
in Ref. 0, is based on a SrTiOs (STO) single crystal 
gate dielectric, the substrate itself. Due to its large low 
temperature dielectric constant, e, and large achievable 
polarizations, STO is a particularly interesting gate in- 
sulator. Field effect devices using a STO thin film gate 
insulator have thus been studied extensively Here, 
thin 100 /im thick, or thinned 500 /im commercial sub- 
strates have been used. A sketch of the thinned device is 
shown in the upper inset of FigQ The superconducting 
NBCO thin films are first grown by off- axis magnetron 
sputtering onto (001) (bare or etched) STO substrates 
heated to about 730°C. During cooling, a typical O2 pres- 
sure of 670 Torr is used to obtain optimally doped films. 
To control the initial doping level of the films, the oxygen 
cooling pressure is modified and lowered down to 5 mTorr 
for the most underdoped films. Gold is sputtered in-situ 
at room temperature to improve the contact resistances 
and the whole structure is protected by an amorphous 
NBCO layer also deposited in-situ. X-ray diffraction al- 
lowed us to determine precisely the film thickness down 
to three unit cell thick films. After deposition, the sam- 
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FIG. 1: Resistance versus temperature for a 3-4 unit cell thick 
NBCO film (sample 1) for different applied gate voltages; the 
right scale is the corresponding resistivity. The lower inset 
shows the resistance as a function of temperature for sample 
5, which is also 3-4 unit cell thick, for different applied gate 
voltages (see text for details). The upper inset shows a sketch 
of the field effect device, based on a thinned STO single crystal 
gate insulator. 

pies are photolithographically patterned using ion milling 
and a gold electrode is deposited on the backside of the 
samples, facing the central part of the superconducting 
path. The measured path is 600 (length) x 500 (width) 
fim^. The patterning process does not affect Tc- All the 
thin films presented here are 3-4 unit cell thick. 

Figfflshows resistance versus temperature for sample 1 
in the tail of the transitions, for different voltages applied 
across the gate dielectric: 0, -20, -50, -100 and -200V. 
Negative voltages correspond to a negative potential ap- 
plied to the gate and a positive one to the superconduct- 
ing path, resulting in hole doping of the oxide channel 
and, as expected, in a "shift" to higher temperatures 
of the resistive transition. Resistance measurements are 
performed using a standard four point technique while 
a voltage, Vg, is applied to the gate. During measure- 
ments, gate leakage currents were kept below a few nA, 
while the current used for resistance measurements was 
typically between 1 and 10/zA. The lower inset of FigQ] 
shows resistance versus temperature for sample 5 and for 
gate voltages of 400, 200, 100, 50, 20, 10, 0, -5, -20 , 
-50, -100, and -200 V. At 4.2K, a large increase in re- 
sistance is observed while applying positive gate voltages, 
effectively reducing the hole concentration. 

To estimate we explore the consistency of our resis- 
tivity data with the expected behavior for a Kosterlitz- 
Thouless (KT) transition, namely p = po exp{—bt~^^^), 
where t — T /Tc — 1 and po and b are material dependent 
but temperature independent parameters 0- Accord- 
ingly, consistency with the KT-scenario is established 
when the plot (d In i?/dT)~^^^ vs. t exhibits, near i = 0, 
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FIG. 2: (dlni?/dr)"^''^ vs. T for samples 1 and 5 at various 
gate voltages. The solid lines indicate the consistency with 
the linear KT relationship. Tc — Tkt is determined by the 
condition (dlni?/dr)~^''^ = at Tkt. 

linear behavior and Tc — Tkt is determined by the con- 
dition (dlni?/dT)"^/^ = at Tkt- A glance at FigH 
reveals that the characteristic KT-behavior is well con- 
firmed in an intermediate regime. In sample 1, it is the 
noise level which limits the accessible regime, while in 
sample 5 a rounded transition appears to occur, a sig- 
nature that the correlation length divergence is limited 
by the size of the homogenous domains. Indeed, above 
Tkt the correlation length increases exponentially as 
^ = ^0 exp(6i^^/^). An estimate of the size of the homo- 
geneous regions for sample 5 for T ~ 7.5K, Tkt — 5.67K 
and 5 = 2.63 (values appropriate for sample 5 a,t V — 0), 
leads to ^ (7.5K) /^o — 100. Remarkably enough, with 
^0 ranging from 10 to lOOA, the size of the homogeneous 
domains exceeds 1000 A. 

The consistency with a KT-transition discussed above 
enables us to estimate Tc ~ Tkt and determine its gate- 
voltage dependence. The resulting Tkt estimates for var- 
ious samples versus gate voltage V are plotted in Fig|31i. 
As can be seen, Tkt (V) is nonlinear and, as indicated by 

1 /2 

the respective solid curves, Tkt (V) ~ Tkt (0)-f a |y| 

To correlate the field-induced Tkt modulations to the 
field-induced areal carrier density, n2D , we measured the 
accumulated charge Q. Since the dielectric constant of 
STO depends on temperature and applied electric field, 
we measured the field and temperature dependence of 
the capacitance The induced areal charge density, 

a = eAn2D, at a given temperature was measured in 
terms of the capacitance of the device as a function of 
the applied gate voltage, using an LCR meter (Agilent 
4284A) with an auto-balancing bridge method, and/or by 
measuring the charge fiow during loading using an elec- 
trometer (Keithley 6514). By ramping the voltage across 
the dielectric, the LCR-meter measures the "local" ca- 
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FIG. 3: (a) Right scale, Tkt versus gate voltage V . The 
solid lines are Tkt (V) ~ Tkt (0) + a \V\^^^ with a = 0.54, 
0.62, 0.32, 0.34 for samples 1, 2, 5 and 6. Left scale, dQ/dV 
versus V for a lOO/im thick STO measured at 4.2 K using an 
electrometer (open circles), C{V) measured using an LCR- 
meter (open squares), (h) ATkt vs. An2D for all the samples. 
The solid line is Eq.Q. 

pacitance C{V) at a given voltage while the electrome- 
ter measures dQ for a given change in voltage j22i|. The 
two measurements agree quantitatively as can be seen on 
FigEli where dQ/dV and C{V) {dQ = C(y)dV), mea- 
sured at 4.2K, are plotted versus V for a 100/xm thick 
STO substrate with two parallel 20 mm^ square gold elec- 
trodes. The corresponding field induced charge density 
(j{V) at a given voltage and temperature is then obtained 
from (7{V) = 1/5 C{V)dV, where S is the surface of 
the gate contact. This quantitative estimate of a allows 
a change in Tkt to be related to a change in the areal 
carrier density — cr/e. The result shown in Fig|3|3 
reveals the intrinsic linear relationship 

ATkt = 1-3 IQ-^AnzD, (1) 

where ATkt = Tkt{V) - Tkt{0) IH and Ariao is ex- 
pressed in cm~^. This is a novel and central result of our 
paper. 

Together with the quantum counterpart of the Nelson- 
Kosterlitz relation Tkt cx A~j^(Tr-t), it imphes that 
Tkt , Aob(O) and n2D are universally related by 

Tkt oc n2D oc -2^— , (2) 

with non-universal factors of proportionality. This rela- 
tionship also confirms the theoretical predictions for a 
2D-QSI transition. Indeed, the scahng theory of quan- 
tum critical phenomena 0, 0| predicts that close to a 
QSI transition Tc scales as n|{^, where z is the dynamic 
and u the critical exponent of the zero temperature in- 
plane correlation length, ^qj, oc n^p. Thus, Eq.(|2Jl re- 
veals that 2j7 ~ 1, the signature of a QSI transition in 



FIG. 4: Uo/ks vs. H for different applied voltages for sample 
1. A linear relation between Uo and log(_ff) is observed with 
a slope which depends on the applied electric field. Inset: 
log(p) as a function of 1/T for sample 2 and for three differ- 
ent electrical fields at O.IT. The zero temperature activation 
energy Uo at a given applied electrical field is obtained from 
the slope of the Arrhenius plot. 

2D [mii. Since close to a QSI transition Aq;,(0) scales 

as 1/A ^^(0 ) oc n^^'^'^ [D is the system dimensional- 
ity) it follows that ATkt oc An2D does not only 
uncover a 2D-QSI transition with zF ~ 1, but also im- 
plies that Tkt , Aq;,(0) and are universally related 
by Eq.(|2Jl. However, this relationship differs drastically 
from Tc oc 1/Aa,b,c(0), derived from penetration depth 
measurements on YBa2Cu3 06+a; single crystals 0,01 and 
films , where Tc was reduced by chemical substitution. 

To substantiate our main result further, we explore 
the temperature, electric and magnetic field dependence 
of the resistive transition, extracting the activation en- 
ergy, U, for vortex motion in the liquid vortex phase. 
In the vortex 2D-limit the resistance in a field is ac- 
tivated with activation energy proportional to — In(iJ). 
The measurements have been performed with the mag- 
netic field applied perpendicular to the a6-plane, ramping 
the temperature slowly and measuring the sample resis- 
tance at a given magnetic and electric field. Inset of 
Fig0| shows Arrhenius plots of the resistivity, log(p) as a 
function of l/T, for sample 2 and for three different elec- 
trical fields at H — O.IT. The observed linear relation 
between log(p) and 1/T singles out thermally activated 
flux flow, where p{H, T) = p„ exp(-J7(iJ, T)JkBT), with 
U{H,T) = 2Uo{H) (1 - T/Tc) for T ~ Tc |lj As can 
be seen, U becomes larger for electrical fields raising the 
number of holes and Tc. We note that an electrostatic 
modulation of the activation energies was also observed 
in reference From such Arrhenius plots we estimated 
Uo{H)/kB for sample 1 at different apphed voltages and 
different magnetic fields. As can be seen in Fig01 at 
a given voltage, we observe between 0.1 and 7 T the 
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characteristic 2D logarithmic field dependence Uo{H) = 
—a ln(_ff ) + /3, in agreement with previous measurements 
on thin films, superlattices, and the hi ghly anisotropic 
kappa-(BEDT-TTF)2Cu(NCS)2 [11 IM llllJil il 13 . 
Note that in 3D J7 oc i/^^/^ is expectedjig"^ Further- 
more we find that Tc — Tkt is proportional to Uo{H)/kB 
for every value of the magnetic field. Since U is propor- 
tional to l/A^^ in D = 2 (and £> = 3) jll El, this 
proportionality is also consistent with our main result 

In summary, we explored the relationships between 
Tc, the mobile areal carrier density, n2D, and the zero 
temperature in-plane penetration depth, Aah(O), for very 
thin underdoped NBCO films near the superconductor 
to insulator transition by means of the electric field ef- 
fect technique. Together with the observed behavior of 
the resistive transition we established remarkable consis- 
tency with 2D critical behavior and a quantum supercon- 
ductor to insulator transition, characterized by a linear 
relationship between Tc, n2D and 1/A^j(0). This result 
also implies that isotope or pressure effects on Tc, n2D 
and Aafc(O) are related accordingly. 
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